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CARBON MONOXIDE OXIDATION RATES COMPUTED FOR AUTOMOBILE 


EXHAUST MANIFOLD REACTOR CONDITIONS 
by Richard S. Brokaw and David A. Bittker 
Lewis Research Center 

SUMMARY 

Carbon monoxide oxidation rates have been computed by integrating differential equa- 
tions for a system of 29 reversible chemical reactions. At temperatures above about 
1030 K (1400° F) oxidation is rapid and is complete in times short in comparison to the 
residence times in exhaust manifold reactors. This suggests that such reactors are mix- 
ing rate limited, not reaction rate limited. At temperatures below 1030 K rates may be 
considerably slower. At these lower temperatures the reaction is catalyzed by nitric 
oxide. Thus, elimination of nitric oxide from automotive exhausts might degrade the low 
temperature performance of exhaust manifold reactors. 


INTRODUCTION 

The thermal reactor is a candidate system for reducing the carbon monoxide and 
hydrocarbon content of automobile engine exhausts. In this scheme the exhaust manifold 
is replaced by a large, well -insulated volume. Secondary air is injected at the cylinder 
exhaust ports and reacts with the pollutants in this thermal reactor. Such reactors are 
effective in reducing carbon monoxide and hydrocarbons in the exhaust; however, ques- 
tions as to cost and durability remain (ref. 1). 

The thermal reactor provides volume both for mixing of the exhaust gas with added 
air and also for subsequent reaction. It is appropriate to inquire as to what extent the 
cleanup of pollutants is limited by mixing rate and to what extent chemical reaction rates 
are limiting. With the rate limiting process identified, a more rational design of the 
thermal reactor should be possible. 

Carbon monoxide and hydrogen are the major combustibles in the untreated exhaust 
gas and may amount to several percent by volume. (Hydrogen is of no concern as a pol- 
lutant.) Hydrocarbon concentrations are considerably smaller, amounting to several 


hundredths of a percent, or several hundred parts per million. While the kinetics of 
hydrocarbon oxidation is not known in detail, the kinetics of the oxidation of mixtures of 
carbon monoxide and hydrogen is well established. Thus, calculations of the rate of 
carbon monoxide oxidation can be made for conditions appropriate to the automotive ther- 
mal reactor. Since the concentration of hydrocarbons is small, their presence should 
have little or no effect on the overall behavior of the combustion process. 

Conditions which are conducive to rapid carbon monoxide and hydrogen oxidation - 
namely large concentrations of hydrogen atoms, oxygen atoms, and hydroxyl radicals - 
will also bring about rapid oxidation of the hydrocarbons. Thus, the residence times 
that are estimated for the cleanup of carbon monoxide should also apply for the cleanup 
of hydrocarbons in exhaust manifold reactors. 

This report presents calculated carbon monoxide oxidation rates for conditions be- 
lieved appropriate to the automotive thermal reactor. The calculations suggest that tem- 
peratures in excess of 1030 K, or about 1400° F, are required for rapid cleanup of pol- 
lutants and that present reactors may be mixing rather than reaction rate limited. 


COMPUTATIONAL DETAILS 
Chemical Kinetics 

Calculations were carried out assuming a system of 29 reversible chemical reactions 
involving compounds of hydrogen, oxygen, carbon, and nitrogen. In general, rate con- 
stants were taken from as few sources as possible in hopes that some internal self- 
consistency might be maintained. A critical review of the literature was not attempted; 
rather, we have tried to draw on the reviews of others insofar as possible. 

We list in appendix A the individual reactions included in our scheme and indicate the 
values and sources of the assumed rate constants, together with comments on the impor- 
tance and role of the particular step. 


Numerical Computation 

Concentrations of the various species as a function of time were determined by sum- 
ming the extent of reaction of a large number of simultaneous and competing chemical 
reactions. In these, the individual rates are dependent on concentration and temperature 
and the temperature, in turn, depends on changes in composition. This requires the 
integration, by machine, of a set of first order, nonlinear differential equations. The 
system was assumed to be adiabatic and to obey the ideal gas law. 
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The solution is part of a general chemical kinetics computer program which solves 
many types of complex chemical reaction problems in flowing and static systems. The 
integration technique is particularly suited to handle efficiently the extremes of both slow 
and fast chemical reaction. A brief discussion of the computer program is given in ref- 
erence 2. 


Input Conditions 

The assumed exhaust gas compositions were appropriate for a typical hydrocarbon 
fuel with a H/C ratio of 1. 9 burned in an engine to two air -fuel ratios - A/F = 12, typ- 
ical of an idling engine, and A/F = 14, a cruise condition. These compositions were 
estimated from Chart 8 of reference 3. Calculations were carried out for these gas mix- 
tures diluted with additional air (21 percent 0 2 > 79 percent N 2 ) first to stoichiometric, 
A/F = 14.64, and then to a substantially leaner ratio, A/F = 17. Thus, four gas mixtures 
were considered. The undiluted and diluted compositions are presented in table I. 

For each mixture kinetic calculations were run at 100° intervals from 900 to 1300 K 
(1160° to 1880° F); in all cases the pressure in the exhaust manifold reactor was as- 
sumed to be atmospheric. Calculations were generally extended to 20 milliseconds which 
is typical of the shorter residence times for the thermal reactor. 


RESULTS AND DISCUSSION 

The computations contain a great deal of detailed information as to the variation of 

gas composition and temperature with time. This is illustrated in figures 1(a) and (b), 

which are semilogarithmic plots of mole fraction as a function of time. This particular 

calculation is for an idle condition (initial A/F = 12) diluted to stoichiometric (A/F = 

14. 64) with an initial temperature of 900 K. In this case, and in all but one of the others, 
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an initial concentration of 1000 ppm of NO was assumed (10 ” mole fraction). The loga- 
rithm of temperature rise is shown as a dashed curve in figure 1(a). 

Initially small concentrations of atoms and free radicals are produced, principally 
oxygen atoms by means of reaction 17 (see appendix A). Chain branching then sets in by 
way of reactions such as 2, 3, and 9 which involve an atom or free radical and a stable 
molecule: 


+ 0 2 - OH + O 

(2) 

+ H 2 ► OH + H 

(3) 
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20H 


(9) 


o + h 2 o 

The chain may be propagated and the identities of the chain carriers can be changed by 
reactions such as 


OH + H 2 

H 2 0 + H 

( 1 ) 

H + 0 2 + M — 

- H0 2 + M 

(4) 

OH + CO C0 2 + H 

(5) 

H°2 + « 2 

H 2 0 2 + H 

(ID 

H0 2 + NO -N0 2 + OH 

( 21 ) 

H + NO + M — 

— HNO + M 

(24) 

These reactions cause an exponential rise in the concentrations of H, OH, NO 2 , and HNO ; 

as is apparent in figure 1(a). (The H0 2 and O concentrations also rise exponentially, but 

have not been plotted to prevent fig. 1 (a) from becoming unduly cluttered.) 

-5 

The mole fractions of N0 2 and H0 2 rise to about 10 and then decline; the mole 
fractions of HNO and H 2 0 2 only rise to 10"®. These species are limited to these small 
concentrations because they are destroyed by rapid bimolecular processes such as 

OH + H0 2 * 

H 2 ° + °2 

(7) 

0 + no 2 * 

NO + 0 2 

(23) 

OH + HNO - 

H 2 0 + NO 

(26) 

OH + H 2 0 2 * 

h 2 o + ho 2 

(16) 

3 

On the other hand, the mole fractions of H, OH, and O rise as high as 10 
species are eliminated by slower termolecular reactions: 

because these 

H + H + M 

~ H 2 + M 

(-18) 
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H + OH + M 


(19) 


H 2 0 + M 

0 + 0 + M 0 2 + M (20) 

The rapid cleanup of carbon monoxide occurs by means of 

CO + OH -C0 2 + H (5) 

_ 4 

When the hydroxyl radical mole fraction exceeds 10 . 

All of the computations included detailed concentration histories such as figures 1(a) 
and (b). However, since our interest here is in carbon monoxide oxidation only, plots of 
carbon monoxide concentration as a function of time will be presented in figures 2(a) to 
5(a) as linear plots and figures 2(b) to 5(b) as semilogarthmic plots. 

In general, there is an induction period, followed by a period of more or less rapid 
consumption of carbon monoxide. The induction period is decreased as the initial tem- 
perature increases. This is a period of rapid buildup of the concentrations of H, OH, 
and O. In an actual exhaust manifold reactor it seems likely that the exhaust gases from 
the cylinders would have large concentrations of atoms and free radicals left over from 
combustion in the cylinders, which would eliminate the induction period. Indeed, it is 
common practice to add air as close to the exhaust valves as possible; perhaps this is 
because high concentrations of atoms and free radicals effectively initiate the oxidation 
which is later completed in the reactor. 

Following the induction period there is more or less rapid oxidation. For all the 
calculations at 1100, 1200, and 1300 K (1520°, 1700°, and 1880° F) this oxidation is in- 
variably rapid, with carbon monoxide concentrations halved in a millisecond or less. 
Reaction rates are higher in the mixtures corresponding to the idle condition (figs. 2(a), 
2(b), 3(a), and 3(b)), due to the higher initial carbon monoxide and hydrogen concentra- 
tions. The final cleanup of carbon monoxide is faster when there is excess air - compare 
figure 2(b) with 3(b) and 4(b) with 5(b). (The late stages in the oxidation are seen more 
clearly in the semilogarithmic plots. ) 

At temperatures of 900 and 1000 K (1160° and 1340° F) reaction may be very much 
slower - see especially figures 4(a), 4(b), 5(a), and 5(b). At these temperatures the rate 
of the chain breaking reaction 


H + 0 2 + M H0 2 + M 

is more than twice the rate of chain branching reaction 


(4) 


5 



OH + O 


( 2 ) 


h + o 2 

that is, k^[M] > 2kg. Under these circumstances hydrogen atoms tend to react to form 
the relatively unreactive hydroperoxyl radicals, H0 2 , rather than the much more re- 
active oxygen atoms and hydroxyl radicals. Indeed, if there were no nitric oxide pre- 
sent, reaction would be very slow - indetectable within the time scale of these calcula- 
tions. Thus, we would find ignition and reaction at temperatures where 2kg ^ k^[M] and 
negligible reaction where k^[M] > 2k 2 ; in other words, there is a limit temperature - a 
sort of chemical switch. For these mixtures this limit temperature ranges from 1016 to 
1037 K, roughly 1400° F. Limit temperatures are included in table I. 

However, it is well known that traces of nitric oxide sensitize the ignition of 
hydrogen -oxygen mixtures at temperatures below the limit. This is due to the reaction 

H0 2 + NO -N0 2 + OH (21) 

which transforms the relatively inert perhydroxyl radical into a reactive hydroxyl radi- 
cal. Nitric oxide is regenerated, largely through the reaction 

N0 2 + H NO + OH (22) 

The temperature history of one such low temperature mixture (initial A/F = 12, 
diluted to stoichiometric, initial temperature = 900 K) is shown as a dashed curve in fig- 
ure 1(a). The limit condition, 2kg = k^[M], is achieved at 1025 K after about 3. 1 milli- 
seconds. At that time the H0 2 concentration maximizes (fig. 1(b)) and the hydrogen per- 
oxide concentration experiences a local minimum. Rapid carbon monoxide oxidation fol- 
lows immediately. 

The effect of nitric oxide concentration on the ignition behavior is shown as a dashed 
curve in figure 2(a), a calculation with 0. 05 percent (500 ppm) NO, rather than the 
0. 1 percent (1000 ppm) assumed for the other calculations. The principal effect of re- 
ducing the NO concentration seems to be to increase the induction period by about 1 milli- 
second. A similar change might be expected if the rate constant for reaction (21) (which 
is a guess) were halved. But the nitric oxide concentration should have little effect on 
the course of the calculations at 1100, 1200, and 1300 K as these are above the 1025 K 
limit for this mixture. 

It is difficult to make a definitive statement as to the extent uncertainties in other 
rate constants introduce uncertainty into the computed carbon monoxide oxidation rates. 
However, the computer printout indicates that under many circumstances reactions (1), 
(2), (4), and (5) are important. These rate constants are well known, with uncertainties 
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of 40 percent or less, and one would expect the errors in computed oxidation rates to be 
no larger. 


CONCLUSIONS 

These calculations of carbon monoxide oxidation rates for automobile exhaust mani- 
fold reactor conditions lead to the following conclusions: 

1. At temperatures in excess of about 1030 K (1400° F) oxidation is rapid and pro- 
ceeds essentially to completion in 5 milliseconds or less. This is short in comparison 
to the residence times of 20 milliseconds or more which are typical for exhaust manifold 
reactors, and suggests that such reactors are limited by mixing rates, not reaction rates. 

2. At temperatures below ~1030 K (1400° F) reaction rates may be considerably 
slower. The reaction is catalyzed by nitric oxide and would be extremely slow in the ab- 
sence of this substance. Thus reduction of nitric oxide to proposed levels of about 

200 ppm might degrade the low temperature performance of exhaust manifold reactors. 
Low temperature oxidation rates depend crucially on the rate of the reaction 

HOg + NO -OH + NOg; this rate has been estimated in the absence of experimental 

data. Thus, the low temperature calculations probably indicate only qualitative trends. 

3. Calculations were carried out for exhaust gas compositions corresponding to idle 
and cruise conditions, with dilution air added to bring the compositions to stoichiometric, 
and to a leaner fuel-air ratio (A/F = 17). The mixtures diluted to stoichiometric had 
residual carbon monoxide concentrations of 0.2 to 0. 5 percent, when the initial temper- 
ature was 1100 K or above. On the other hand, mixtures diluted to A/F = 17 had carbon 
monoxide concentrations of 0. 05 percent in 5 milliseconds or less when initial tempera- 
tures were 1100 K or higher. Thus excess air definitely promotes the final cleanup of 
the last traces of carbon monoxide. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 21, 1970, 

129-01. 


7 



APPENDIX A - SPECIFIC REACTION RATE CONSTANTS 


This appendix lists individual reactions included in the calculations and indicates 

values and sources of the assumed rate constants. Bimolecular rate constants are ex- 

3 11 fi 2 

pressed in units of cm mole - sec" and termolecular rate constants are in cm mole' 

sec - *. Temperatures are in degrees Kelvin and activation energies are in calories. 

1. OH + H 2 -H 2 0 + H kj = 2. 3X10 13 exp(-5200/RT) 

This is an important chain propagation step in the branched -chain hydrogen -oxygen 
ignition scheme. The rate is very well established; the rate constant is taken from ref- 
erence 4. 

2. H + 0 2 OH + O k 2 = 2. 04X10 4 exp(-16 500/RT) 

This rate constant, which is from reference 4, is the crucial branch in the hydrogen 
ignition scheme. A recent shock tube study suggests rate constants may be 35 percent 
smaller in the 1200 to 1400 K range (ref. 5). 

3. O + H 2 OH + H k 3 = 4xl0 13 exp(-10 200/RT) 

This is another branch in the H 2 -0 2 reaction. This rate constant, from reference 4, 
may be 15 to 25 percent too large (refs. 5 and 6). 

4. H + 0 2 + M - H0 2 + M k 4 = 10 15 exp(+1300/RT) 

The rate constant shown is for M = argon. For other third bodies the rate constant 
was multiplied by the following factors: H 2 0 , 30; H 2 , 5; 0 2 , 1.75; CO, 2; C0 2 , 7.35; 
N 2 , 2. 15; and all others, 1. This important chain termination gives rise to the second 
limit for the H 2 -0 2 reaction. The rate constant is from reference 4. The relative third 
body efficiencies are from Lewis and von Elbe (ref. 7), except for H 2 0 (ref. 8) and CO 
(a guess). 

5. CO + OH -C0 2 + H k 5 = 6. 6x10 11 exp(-1030/RT) 

This is the principal path for the conversion of carbon monoxide to carbon dioxide. 
The rate constant is taken from Schofield (ref. 9). Recent shock tube results suggest a 
value lower by perhaps 35 percent in the range 1300 to 1900 K (ref. 5). 

6. H + H0 2 20H kg = 7xl0 13 

This rate constant is from reference 4. It is close to collision frequency so that the 
activation energy must be very small. 

7. OH + H0 2 ~H 2 0 + 0 2 k 7 = 6xl0 12 

Kaufman (ref. 10) suggests that the rate constant is this large or larger. 
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8. O + H0 2 -OH + 0 2 kg = 6x10 12 

Again, a lower limit proposed by Kaufman. 

9.0 + H 2 0 20H k g = 8. 4xl0 13 exp(-18 OOO/RT) 

Rate constant suggested by Schofield (ref. 9). 

10. H + H0 2 H 2 + 0 2 k 1Q = 2. 3xl0 13 

The reverse of this reaction provides a source of free radicals from H 2 and Og. 
Clyne and Thrush (ref. 11) suggest that kg/k^g ~ 3; this combined with Baldwin’s rate 
(ref. 4) for kg leads to k^g ~2.3xl0* 3 . 

11. H0 2 + H 2 H 2 0 2 + H k n = 1. 66xl0 13 exp(-25 000/RT) 

This is a slow reaction regenerating hydrogen atoms from the relatively inert H0 2 
radicals. This rate constant is from Baldwin (ref. 4). 

12. H 2 0 2 + M -20H + M k l2 = 3. 19xl0 17 exp(-47 000/RT) 

This rate constant is for M = nitrogen from reference 4. For other third bodies the 
rate constant was multiplied by the following factors : HgO, 6; 0 2 , 0.78; H 2 , 2.3; H 2 0 2 , 
6.6; and all others, 1. These factors are from Baldwin and Brattan (ref. 12) except for 
H 2 (private communication from Professor R. R. Baldwin). 

13. 2H0 2 ~H 2 0 2 + 0 2 k 13 = 1. 8x10 12 

This rate constant is from reference 4. 

14. H + H 2 0 2 H 2 0 + OH k 14 = 4. 16xl0 14 exp(-9000/RT) 

This rate constant is from reference 4. 

15. O + H 2 0 2 OH + H0 2 k 15 = 8X10 13 exp(-1000/RT) 

This rate constant is estimated by analogy with reaction 9. The activation energy of 
reaction 9 exceeds the endothermicity by 1000 calories. Since reaction 15 is exothermic, 
we assume an activation energy of 1000 calories with a preexponential constant similar to 
that for reaction 9. 

16. OH + H 2 0 2 H 2 0 + H0 2 k lg = 3. 6X10 12 

This rate constant is from reference 4. This is in reasonable agreement with more 
recent results (ref. 13). 

17. CO + 0 2 -C0 2 + O k 17 = 2. 5X10 12 exp(-48 000/RT) 

This reaction between carbon monoxide and oxygen serves as an initiating source 
of atomic oxygen. The rate constants shown, from reference 14, may be too low, by an 
order of magnitude or more (ref. 5). 
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18. H 2 + M 2H + M k 18 = 1. 12xl0 13 T 1//2 exp(-92 600/RT) 

The rate constant shown is for M = H 2 (ref. 15). This value was used on the notion 
that several of the components of exhaust gas mixtures - H 2 , HgO, C0 2 , etc. - should 
be considerably more effective than argon in dissociating hydrogen molecules and re- 
combining hydrogen atoms. A perusal of recent information on relative efficiencies 
(ref. 16) suggests a rate lower by 40 percent might be more realistic for typical auto- 
motive exhaust gas compositions. 

19. H + OH + M H 2 0 + M k 19 = 10 19 /T 

This rate is for M = Ar. The rate constant was estimated in the following manner: 

The rate of this recombination is about ten times as fast as the rate of hydrogen atom 

recombination (refs. 16 and 17), the reverse of reaction 18. And a shock tube study 
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(ref. 18) gives k 1fi = 10 /T for M = Ar. More recent work (ref. 19) gives k 1Q = 
o o p c “ 1 o iy 

7. 5x10^ /T for M = argon over a wide temperature range, which gives rate constants 
within 40 percent of those we have used in the temperature range of interest here (900 to 
1300 K). However, since nitrogen is 1.6 times as effective as argon and water is twenty 
times as effective (ref. 16), the rate of this reaction might be roughly four times faster 
in exhaust gas mixtures. 

20. 20 + M ►Og+M k 2Q = 8. 1 Sx10 18 /T 1,22 

This rate, from reference 20, is for M = oxygen. A recent critical survey (ref. 21) 
suggests rates 30 to 40 percent lower in the range 900 to 1300 K. 

21. NO + H0 2 ► N0 2 + OH k 21 = 10 13 

This reaction converts the relatively unreactive hydroperoxyl radical to the reactive 
hydroxyl radical. It is responsible for the NO sensitized ignition at temperatures below 
the second limit, conditions where k^[M] > 2kg. We have assumed a moderately fast 
bimolecular rate, similar to reactions 6, 7, 8, and 10. 

22. N0 2 + H ►NO + OH k 22 = 7.2X10 14 exp(-1930/RT) 

This rate constant is from reference 9. 

23. O + N0 2 ►NO + 0 2 k 23 = 1. 9xl0 13 exp(-1060/RT) 

This rate constant is from reference 9. 

24. H + NO + M ►HNO + M k 24 = 4xl0 15 exp(+600/RT) 

The rate constant shown is for M = argon. It was obtained by combining the 293 K 
datum of Hartley and Thrush (ref. 22) with the temperature dependence suggested by 
Schofield (ref. 9). For other third bodies the rate constant was multiplied by the follow- 
ingfactors: H 2 0, 6.1; H 2 , 1.87; C0 2 , 2.05; all others, 1. These factors are from 
reference 22. 
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25. H + HNO ^H 2 + NO k 25 = 5xl0 13 

The rate constant used was inadvertently taken a factor of ten larger than the value 
reported in reference 23. A repeat calculation was run for a mixture at 900 K, where 
this reaction should be most important. There was negligible change in the concentration 
profiles except that the HNO concentrations were increased by as much as 40 percent. 

But at most, the HNO concentration barely exceeded 1 ppm. 

26. OH + HNO -H 2 0 + NO k^g = 3.6X10 13 

Rate constant is from reference 23. 

27. O + HNO ►OH + NO k 2? = 3xl0 13 

This estimate obtained by analogy with reaction 26. 

28. H0 2 + NO ►HNO + 0 2 k 2g = 10 13 

This rate is assumed the same as reaction 21. 

29. O + NO + M ►NOg + M k 29 = 9.4xl0 14 exp(+1930/RT) 

The rate shown for M = argon is from reference 9. For other third bodies the 
rate shown was multiplied by the following factors: H 2 0, 6.3; 0 2 , 1.0; N 2 , 1.55; C0 2 , 
2.17; and all others, 1. 

Rate constants for the reverse of these 29 reactions were obtained by invoking de- 
tailed balancing; in other words, the rate constant for each reverse reaction was taken 
as the quotient of the forward rate constant and the equilibrium constant. Thermody- 
namic data were taken from the JANAF thermochemical tables. 


11 



REFERENCES 


1. Cantwell, E. N. ; Roselund, I. T. ; Barth, W. J. ; and Ross, S. W. : Recent De- 

velopments in Exhaust Manifold Reactor Systems. Presented at SAE meeting, 
Detroit, Mich. , Jan. 13-15, 1970. 

2. Bittker, David A. : Chemical Kinetic Computations in Multireaction Systems. Pro- 

ceedings of NASA Conference on Kinetics and Thermodynamics in High Temper- 
ature Gases. NASASP-239, 1970. 

3. Eltinge, Lamont: Fuel -Air Ratio Distribution from Exhaust Gas Composition. Fuels 

and Lubricants. SAE, 1968, pp. 625-653. 

4. Baldwin, R. R. ; Jackson, D. ; Walker, R. W. ; and Webster, S. J. : Interpretation 

of the Slow Reaction and Second Limit of Hydrogen Oxygen Mixtures by Computer 
Methods. Trans. Faraday So c. , vol, 63, 1967, pp. 1676-1686, especially table 12, 
p. 1685. 

5. Brabbs, T. A.; Belles, F. E. ; and Brokaw, R. S.: Shock Tube Measurements of 

Specific Reaction Rates in the Branched Chain H 2 -CO-O 2 System. Paper presented 
at 13th Symposium (International) on Combustion, Salt Lake City, Utah, Aug. 23-29, 
1970. 

6. Westenberg, A. A. ; and de Haas, N. : Reinvestigation of the Rate Coefficients for 

O + H 9 and O + CH 4 . J. Chem. Phys. , vol. 50, no. 6, Mar. 15, 1969, pp. 2512- 
2516. 

7. Lewis, Bernard; and Von Elbe, Guentner: Combustion, Flames, and Explosions of 

Gases. Academic Press, Inc. , 1951, p. 34. 

8. Getzinger, R. W. ; and Schott, G. L. : Kinetic Studies of Hydroxyl Radicals in Shock 

Waves. V. Recombination via the H + O 2 + M — HC >2 + M Reaction in Lean 
Hydrogen -Oxygen Mixtures. J. Chem. Phys. , vol. 43, no. 9, Nov. 1, 1965, 
pp. 3237-3247. 

9. Schofield, K. : An Evaluation of Kinetic Rate Data for Reactions of Neutrals of At- 

mospheric Interest. Planet. Space Sci. , vol. 15, 1967, pp. 643-670. 

10. Kaufman, Frederick: Aeronomic Reactions Involving Hydrogen. A Review of Recent 

Laboratory Studies. Annales de Geophysique, vol. 20, no. 1, Jan. -Mar. 1964, 
pp. 106-114. 

11. Clyne, M. A. A. ; and Thrush, B. A. : Rates of Elementary Processes in the Chain 

Reaction Between Hydrogen and Oxygen. II. Kinetics of the Reaction of Hydrogen 
Atoms with Molecular Oxygen. Proc. Roy. Soc. (London), Ser. A, vol. 275, 

1963, pp. 559-574. 


12 



12. Baldwin, R. R. ; and B rattan, D. : Homogeneous Gas -Phase Decompositions of Hy- 

drogen Peroxide. Eighth Symposium (International) on Combustion. Williams and 
Wilkins Co., 1962, pp. 110-119. 

13. Baulch, D. L. ; Drysdale, D. D. ; and Lloyd, A. C. : High Temperature Reaction 

Rate Data. Rep. no. 3, Dep. Phys. Chem. , Leeds University, England, Apr. 

1969, pp. 40-43. 

14. Brokaw, R. S. : Ignition Kinetics of the Carbon Monoxide - Oxygen Reaction. 

Eleventh Symposium (International) on Combustion. Combustion Institute, 1967, 
pp. 1063-1073. 

15. Myerson, A. L. ; and Watt, W. S. : Atom -Formation Rates behind Shock Waves in 

Hydrogen and the Effect of Added Oxygen. J. Chem. Phys., vol. 49, no. 1, 

July 1, 1968, pp. 425-433. 

16. Getzinger, R. W. ; and Blair, L. S. : Recombination in the Hydrogen-Oxygen 

Reaction: A Shock Tube Study with Nitrogen and Water Vapour as Third Bodies. 
Combustion and Flame, vol. 13, no. 3, June 1969, pp. 271-284. 

17. Schott, Garry L. ; and Bird, Paul F. : Kinetic Studies of Hydroxyl Radicals in Shock 

Waves. IV. Recombination Rates in Rich Hydrogen -Oxygen Mixtures. J. Chem. 
Phys., vol. 41, no. 9, Nov. 1, 1964, pp. 2869-2876. 

18. Jacobs, T. A. ; Giedt, R. R. ; and Cohen, Norman: Kinetics of Hydrogen Halides in 

ShockWaves. II. A New Measurement of the Hydrogen Dissociation Rate. J. 

Chem. Phys. , vol. 47, no. 1, July 1, 1967, pp. 54-57. 

19. Homer, J. B. ; and Hurle, I. R. : The Dissociation of Water Vapour behind Shock 

Waves. Proc. Roy. Soc. (London), Ser. A, vol. 314, 1970, pp. 585-598. 

20. Kiefer, John H. ; and Lutz, Robert W. : Recombination of Oxygen Atoms at High 

Temperatures as Measured by Shock -Tube Densitometry. J. Chem. Phys., 
vol. 42, no. 5, Mar. 1, 1965, pp. 1709-1714. 

21. Johnston, Harold S.: Gas Phase Reaction Kinetics of Neutral Oxygen Species. Rep. 

NSRDS-NBS-20, National Bureau of Standards , Sept. 1968. 

22. Hartley, D. B. ; and Thrush, B. A. : The Rate of Elementary Processes in the Chain 

Reaction between Hydrogen and Oxygen. HI. Kinetics of the Combination of Hydro- 
gen Atoms with Nitric Oxide. Proc. Roy. Soc. (London), Ser. A, vol. 297, 1967, 
pp. 520-533. 

23. Halstead, C. J. ; and Jenkins, D. R. : Catalysis of Recombination Reactions in 

Flames by Nitric Oxide. Chem. Phys. Letters, vol. 2, no. 5, Sept. 1968, 

pp. 281-282. 


13 



TABLE I. - ASSUMED EXHAUST GAS COMPOSITIONS 8 


H 2 H 2 0 CO C0 2 0 2 NO 

Idle condition (engine A/F = 12) 

Undiluted 2.19 12.37 5.87 9.46 0.10 

Diluted to stoichiometric 1.84 10.38 4.93 7.94 3.33 . io 

(overall A/F = 14. 64; 
limit, 1025 K) 

Diluted lean (overall A/F = 1.61 9.06 4.30 6.93 5.57 .10 

17; limit, 1016 K) 



Cruise condition (engine A/F = 14) 


Undiluted 0.41 12.70 1.40 12.40 0.10 

Diluted to stoichiometric .39 12.17 1.34 11.89 0.816 .10 

(overall A/F = 14.64; 
limit, 1037 K) 

Diluted lean (overall A/F = .34 10.59 1.17 10.34 3.44 .10 

17; limit, 1028 K) 

a Percent by volume. 
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